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Description 

VOXEL-DRIVEN SPIRAL RECONSTRUCTION FOR CONE BEAM COMPUTED 
TOMOGRAPHY 

The invention relates to a method for generating images in 
computed tomography with the aid of 3D image reconstruction 
method, wherein in order to scan an examination object with the 
aid of a conical beam emanating from a focus and the aid of a 
planar, preferably multirow, detector for detecting the beam, 
the focus is moved on a spiral focal track about the 
examination object, the detector supplying output data, that 
correspond to the detected radiation, and image voxels from the 
scanned examination object are reconstructed from the possibly 
preprocessed output data and reproduce the attenuation 
coefficients of the respective voxel. 

Beams having cone beam geometry are used in the prior art of 
computed tomography with multirow detectors. Reconstructing an 
image volume requires consideration of the cone beam geometry 
in the 3D image reconstruction, the cone beam geometry 
substantially complicating the reconstruction problem. Two 
different groups of image generation methods are known in 
principle, specifically approximate methods and exact methods: 

Approximate methods : 

Approximate methods are distinguished by a high measure of 
practicability and flexibility. The angle of inclination of the 
measuring beams to the axis of rotation (cone angle) is 
considered in an approximate fashion, for which reason the 
error owing to the approximation grows with the cone angle. 
Starting from a certain number of detector rows, each 
approximate method will cause 
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image artefacts. A distinction is made between 2D and 3D 
methods in the case of the approximate methods. 

ASSR [1], AMPR [2] and SMPR [3] may be named as examples of the 
2D method or 2D rebinning. In this case, synthetic projection 
data are approximated from the cone beam data, the geometry of 
the synthetic projection data being selected such that all the 
synthetic measuring beams lie in a flat plane. The object 
distribution in the flat plane can be reconstructed from the 
synthetic projection data with the aid of the conventional 
algorithms of 2D CT. However, with growing cone angles the 
approximation in the rebinning step leads relatively quickly to 
image artefacts. 

In the case of the 3D method, filtering of the projection data 
and a subsequent 3D back projection are undertaken. The 
recorded geometry of cone beam type is taken into consideration 
here exactly in the back projection. Various feasible 
approximations are applied in the filter step. These 
approximation methods set forth below have in common the fact 
that the filtering consists of a ID ramp filter- that is aligned 
in the direction of the projection of the spiral tangent. The 
individual approximation methods differ in the processing of 
redundant data. 

On the one hand, it is possible to average in the axial 
direction, as is known from [4] . Here, after the filter step 
the filtered projection data are back projected into the 
reconstruction volume such that the measured data are 
accumulated in weighted fashion in the process. The weights 
result from the axial distance of the measuring beams from the 
voxel to be reconstructed. 

On the other hand, a reconstruction of axial slices is known 
from [5] . For this purpose, all the measured data that cut the 
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image slice are used for the reconstruction. Data redundancies 
are considered approximately, neglecting the 
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cone angle. Use may be made for this purpose of conventional 
methods such as, for example, the Parker weighting known from 
[6] . The filter step consists in this case of a weighting of 
the projection data in accordance with the present data 
redundancy, followed by convolution using a ID ramp filter. 
This method requires a relatively large detector surface, since 
the projection of the entire reconstruction slice must be 
included on the detector for each focal position. 

Exact methods : 

A comparative study of the most important exact algorithms is 
to be found in [7]. Exact methods consider the recorded 
geometry of cone beam type in a fashion free from error both in 
the filter step and in the 3D back projection. These methods 
achieve good image results independently of the cone angle that 
occurs. However, they are extremely complicated and very 
inflexible in application. For example, a reduced pitch 
accompanied by the employment of data redundancies can be 
realized only to a limited extent, and there is no possibility 
at all of selecting measured data for cardiac imaging with a 
high time resolution. 

The problem now consists in that the image quality achieved and 
the required computational outlay are oppositely oriented as a 
function of the method used. Algorithms having excellent image 
quality in conjunction with large cone angles, that is to say 
for detectors with many detector rows, are essentially intended 
for future scanner generations. It would be desirable to avoid 
exact methods, since these are complicated and inflexible. 

It is therefore an object of the invention to find in computed 
tomography having multirow detectors and beams with cone beam 
geometry a proposal to solve the generation of images that, on 
the one hand constitutes an optimum compromise between 



PCT/EP2004/000754 
2002P13580WOUS " 



- 4 - 



image quality and computational outlay, and on the other hand 
can also be applied sufficiently flexibly. 

This object is achieved by means of the features of the 
independent patent claims. Advantageous developments of the 
invention are the subject matter of dependent claims. 

The inventors have recognized the following: 

Starting from a reconstruction method such as is described in 
[5] , instead of reconstructing entire image slices it is also 
possible to reconstruct the image voxels separately. It is then 
necessary in this case to process data redundancies by using 
methods that permit every arbitrary projection angular range 
that is greater than the required half revolution n. 

What is novel in the invention is the separate reconstruction 
of individual voxels. The projection data are weighted 
separately for approximate consideration of data redundancies 
for each individual voxel. The result of this is that, by 
contrast with the conventional back projection, instead of a 
filtered measured value being smeared along a line over all the 
voxels cut by the line, the measured value filtered for the 
voxel V is accumulated only on the voxel V. 

The present invention offers the following advantages, good 
image quality is expected even for large cone angles, 
flexibility in pitch up to the theoretical maximum value, use 
of all the data redundancies and the possibility of data 
selection, for example, for cardiac imaging with high time 
resolution . 

In accordance with this basic idea, the inventors propose the 
following method for generating images in computed tomography 
with the aid of a 3D image reconstruction method that has at 
least the following method steps: 
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in order to scan an examination object with the aid of a 
conical beam emanating from a focus and the aid of a 
planar, preferably multirow, detector for detecting the 
beam, the focus is moved on a spiral focal track about the 
examination object, the detector supplying output data, 
that correspond to the detected radiation, and 
image voxels from the scanned examination object are 
reconstructed from the possibly preprocessed output data 
and reproduce the attenuation coefficients of the 
respective voxel, 

each image voxel being reconstructed separately from 
projection data that comprise a projection angular range 
of at least 180°, and 

an approximate weighting taking place for each voxel 
considered in order to normalize the projection data used 
relating to the voxel. 

It is advantageous in this case when - as illustrated in Figure 
4 - in order to reconstruct an image voxel V, use is made of 
all the detector data along a straight line that runs through 
the cone beam projection of the voxel V and is aligned in the 
direction of the projection t of the spiral tangent S t . An 
optimum detector employment is achieved thereby. 

Furthermore, the image data of the detector image can be 
subjected to cosine weighting in order thereby to achieve a 
compensation of the oblique irradiation of the focus onto the 
detector . 

If the scanning operation does not have a sufficiently high 
resolution in order for all the required detector data to be 
obtained directly for a voxel being considered, data not 
directly available can be obtained from the available data by 
interpolation from neighboring detector data (detector pixel 
values) . 
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Furthermore, during the weighting for compensating the data 
redundancy (normalization) two measuring beams (Sa, Sb) can be 
regarded as redundant precisely when it holds for the parallel 
coordinates 0, p that: (9 a = 2k-7i + 0 b and p a = p b ) or 
(0 a = (2k + l)-7i + 0 b and p a = -Pb) / where k represents an 
arbitrary natural number, and 0 represents the projection 
angle, and the projection angle 0 represents the sum of focal 
angle a and fan angle (3, and p represents the distance from the 
z axis. 

The redundantly available data can be multiplied by generalized 
Parker weights for the purpose of normalization. The normalized 
data are processed with the aid of a ramp filter modified by an 
advantageous smoothing window. It is advantageous, furthermore, 
when a distance weighting is performed for the 3D back 
projection into the voxel considered. 

Finally, it is to be pointed out that the method according to 
the invention can also be used for cardiac computed tomography 
by selecting, weighting or sorting in a manner known per se the 
measured data in accordance with the movement phases of an 
examined heart, in order thus to obtain tomograms of the heart 
in specific movement phases. 

In accordance with the basic idea of the invention, the 
inventors also propose a CT unit for scanning an examination 
object that is provided with a detector array that is of planar 
design and has a multiplicity of distributed detector elements 
for detecting the rays of the beam from at least one focus, the 
at least one focus being movable relative to the examination 
object on at least one focal track that runs around the 
examination object and a detector array situated opposite, at 
least means for collecting detector data, filtering and 3D back 
projection being provided, and the means for processing the 
measured data being fashioned in such a way that the method 
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according to the invention as outlined above can be carried 
out. Said functional means are 
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preferably implemented at least partially by programs or 
program modules . 

The invention is described in more detail below with reference 
to a preferred exemplary embodiment of a third-generation 
multirow CT unit suitable for carrying out the method according 
to the invention, and with the aid of the figures in which: 



Figure 1 shows multirow unit in a perspective illustration of 
the scanning unit in a schematic illustration of the 
peripheral equipment ; 

Figure 2 shows a longitudinal section through the multirow CT 
from Figure 1; 

Figure 3 shows illustrations of the recording geometry of a 
spiral CT unit; 

Figure 4 shows a schematic of the collection of measured data 
and of the approximate redundancy of measuring beams; 
and 

Figure 5 shows an illustration of the angular interval A(V) . 



Figures 1 and 2 show a partially perspective illustration of a 
third-generation multirow CT unit suitable for carrying out the 
method according to the invention. The measuring arrangement (= 
Gantry) denoted by 1 has an x-ray source 2 with a beam 
diaphragm 3 in front thereof and near the source. The two- 
dimensionally constructed array with a number of rows and 
columns of detector elements forms the detector system 5 and is 
represented in section in figure 2 with a beam diaphragm 6 in 
front thereof and near the detector. For the purpose of better 
clarity, only eight rows Li to L Q of detector elements 4 are 
illustrated in figure 1. However, the detector system 5 can 
also have another, preferably greater, number of rows without 
departing from the scope of the invention. Another two- 
dimensional arrangement of the detectors is likewise also 
possible . 



PCT/EP2 004/0007 54 
2002P13580WOUS " 



- 8 - 



The X-ray source 2 with the radiation diaphragm 3, on the one 
hand, and the detector system 5 with the radiation diaphragm 6, 
on the other hand, are arranged opposite each other on a rotary 
frame 7 in such a way that a pyramidal X-ray beam which, during 
the operation of the CT unit, emanates from the X-ray source 2 
and is collimated by the settable radiation diaphragm 3 and 
whose edge rays are designated by 8 , strikes the detector 
system 5. In the process, the radiation diaphragm 6 can be set 
to correspond to the cross section, set by means of the 
radiation diaphragm 3, of the X-ray beam in such a way that, in 
accordance with different operating modes, only that area of 
the detector system 5 that is struck directly by the X-ray beam 
is exposed. In figures 1 and 2, only eight rows of detector 
elements 4 are used, the further rows indicated by dots are 
covered by the radiation diaphragm 6 and are therefore 
inactive . 

The X-ray beam has a cone angle <p which is the opening angle of 
the X-ray beam in a plane containing the system axis Z and the 
focus F. The opening angle of the X-ray beam in a plane lying 
at right angles to the system axis Z and containing the focus F 
(fan opening angle) is 2/jLix. 

The rotary frame 7 can be set rotating about a system axis 
designated by Z by means of a drive device 22. The system axis 
Z runs parallel to the z-axis of a three-dimensional 
rectangular coordinate system illustrated in figure 1. 

The columns of the detector system 5 likewise run in the 
direction of the z-axis, while the rows run transversely with 
respect to the system axis Z. 
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In order to be able to bring an examination object, for example 
a patient, into the beam path of the X-ray beam, a bearing 
device 9 is provided, which can be displaced parallel to the 
system axis Z, that is to say in the direction of the z-axis, 
specifically in such a way that there is synchronization 
between the rotational movement of the rotary frame 7 and the 
translational movement of the bearing device, with the effect 
of producing a spiral movement of the focus and of a corotating 
detector relative to the examination object. The ratio between 
translational and rotational speeds can be set by a desired 
value for the feed H of the bearing device 9 being selected per 
rotation of the rotary frame. 

It is therefore possible for a volume of an examination object 
which is located on the bearing device 9, to be examined in the 
course of volume scanning, it being possible for the volume 
scanning to be undertaken in the form of spiral scanning in 
such a way that, with simultaneous rotation of the measuring 
unit 1 and translation of the bearing device 9, a multiplicity 
of projections from various projection directions are recorded 
by means of the measuring unit per revolution of the measuring 
unit 1. During the spiral scanning, the focus F of the X-ray 
source is moved relative to the bearing device 9 on a spiral 
path S about the examination object. 

Measured data read out in parallel from the detector elements 
of each active row of the detector system 5 during the spiral 
scanning, and corresponding to the individual projections in 
cone beam geometry, are subjected to digital /analog conversion 
in a data conditioning unit 10, are serialized and transmitted 
to an image computer 11. 

After the measured data have been preprocessed in a 
preprocessing unit 12 of the image computer 11, the resulting 
data stream passes to a volume image reconstruction unit 13, 
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examination object by means of the computing operation 
described later. Both in the preprocessing unit 12 of the image 
computer 11 and in the volume image reconstruction unit 13, 
programs or program modules are stored by means of which the 
method according to the invention can be carried out in order 
to generate CT images with the aid of a 3D image reconstruction 
method . 

The CT images are assembled from voxels (voxel = volume 
element) assembled like a grid, the voxels being assigned to 
the respective image volume, each voxel being assigned a CT 
number in Hounsf ield Units (HU) , and it being possible for the 
individual voxels to be displayed in accordance with a CT 
number/gray value scale in a gray value corresponding to the 
respective CT number. The image volume can be visualized in 
different ways. In a simple form, arbitrary planes can be 
displayed in an arbitrary orientation as tomograms. However, 
there are also more complex methods that visualize the entire 
volume. The shaded surface display (SSD) and volume rendering 
(VR) may be mentioned here by way of example. The SSD 
calculates for a settable threshold value an isosurface for 
which the image volume has the threshold value. The isosurface 
can be displayed on a screen with the aid of methods of 
computer graphics. In the case of VR, each voxel is assigned 
settable optical properties such as opacity and color in 
accordance with its value. Artificial views of the objects thus 
defined are calculated with the aid of methods of computer 
graphics . 

Images are calculated by a visualization unit 15 from the 
volume data reconstructed by the image reconstruction unit 13 
and are displayed on a display unit 16, for example, a monitor 
connected to the image computer 11. 
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The X-ray source 2, for example an X-ray tube, is supplied by a 
generator unit 17 with the requisite voltages and currents, for 
example the tube voltage U. 
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In order to be able to set the latter to the respectively 
requisite values, the generator unit 17 is assigned a control 
unit 18 with a keyboard 19, which permits the values to be set 
as required. It is to be noted in this regard that in addition 
to the keyboard 19 it is also possible to provide other known 
input devices such as mouse, joystick or else voice input. 

In addition, the operation and control of the CT unit apart 
from this is carried out by means of the control unit 18 and 
the keyboard 19, which is illustrated by the fact that the 
control unit 18 is connected to the image computer 11. 

Amongst other things, the number of the active rows of detector 
elements 4, and therefore the position of the radiation 
diaphragms 3 and 6, can be set, for which purpose the control 
unit 18 is connected to the adjustment units 2 0 and 21 assigned 
to the radiation diaphragms 3 and 6. In addition, the rotation 
time x can be set, which is the time needed by the rotary frame 
7 for a complete revolution and which is illustrated by the 
fact that the drive unit 22 associated with the rotary frame 7 
is connected to the control unit 18. 

The recording geometry of the spiral CT unit from figures 1 and 
2 is illustrated schematically in figures 3 and 4, the detector 
being reproduced only in a planar fashion here for the sake of 
simplicity. The radiation source moves with its focus F along a 
spiral track S of pitch h, which corresponds to the bed feed, 
in a fashion rotating about an elongated examination object P, 
and in the direction of the system axis Z. In this process, the 
radiation intensity of the rays penetrating the examination 
object P is measured and collected on the detector D situated 
opposite. A volume V R to be reconstructed is embedded within 
the examination object P. 

In order to reconstruct a volumetric image from the measured 
data, the latter must be filtered, and the filtered data must 
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subsequently be back-projected in three dimensions in order to 
generate a volumetric image of the examination object, the 
volumetric image then representing absorption values, obtained 
from the data, of the voxels V, belonging to the volume of the 
examination object, for the radiation of the beam. 

Thus, consideration of the reconstruction method shows that the 
cone beams from all focal positions that cut an image voxel V 
and are measured on the detector surface contribute to the 
reconstruction of the voxel V. The focal angular range A(V) 
with useful projection data is therefore a coherent angle 
interval, and is given by 

A(v) = [a 0 (v) ,a!(v) ] . 

Figure 5 illustrates this angular interval A(V) . It is assumed 
in this illustration without limitation of generality that the 
spiral track is traversed by focus and detector in a 
counterclockwise sense. The source position F(Oo) is 
distinguished in that the conical projection of the voxel V 
firstly strikes the detector surface D(a 0 ) corresponding to 
the source position F(ao), specifically the upper edge of the 
detector. In the further course of data acquisition, the 
conical projection of the voxel V migrates over the detector 
surface until it leaves the detector at the lower edge in the 
case of the focal position F(a±) . All the projections in the 
angular interval A(V) between Ob and a± contribute to the 
reconstruction of the image voxel V. 

For example, the associated detector images can now be 
processed as follows: 

The first step is for the detector image to be subjected to a 
cosine weighting known per se, in order to take account of the 
oblique radiation present owing to the cone beam geometry. 
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Subsequently, all the detector data are collected per voxel V 
considered, this being done along a straight line that runs 
through the cone beam projection of the voxel V and is aligned 
in the direction of projection t of the spiral tangent S t . If 
these data are not originally available, the required data are 
obtained from neighboring original data by means of an 
interpolation method . 

This is followed by the approximate weighting for the purpose 
of compensating data redundancies. For this purpose, the 
projection angle 0 = a + P and the parallel coordinate p = R F - 
sin (3 (with R F being the radius of the focal track) can be 
determined from the focal angle a and the fan angle P, the cone 
angle being neglected. Two measuring beams S a and S b are 
considered as redundant in this case when it holds that (0 a = 
2k-7i + 0 b and p a = Pb) or (0 a = (2k + 1) -n + 0 b and p a = -p b ) and k 
represents an arbitrary natural number. If redundant data are 
available, they can be multiplied by the generalized Parker 
weights, as is described in [8] . It is to be noted in this 
regard that the measuring beams defined as approximately 
redundant because of the different inclination to the z-axis 
(cone angle) are different and therefore not redundant. 
However, this approximation of the data redundancy can be 
extended relatively far because of the generally relatively 
small cone angle. 

Ramp filtering with a suitable smoothing window can be used to 
filter the values, and the 3D back projection into the voxel V 
can be performed with the aid of a distance weighting known per 
se . 

It may be pointed out in addition that the method according to 
the invention is suitable both for CT units that are operated 
with simultaneously circulating focus/detector pairs, and for 
CT units that have annularly arranged detectors that are fixed 
in the direction of rotation and in the case of which only the 
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direction. 
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